cholesky_rank1_downdate Interface

public interface cholesky_rank1_downdate

Module Procedures

private pure subroutine cholesky_rank1_downdate_dbl(r, u)

Computes the rank 1 downdate to a Cholesky factored matrix such that . This operation only works if the new matrix is positive definite.

Arguments

Type IntentOptional Attributes Name
real(kind=real64), intent(inout), dimension(:,:) :: r

On input, the N-by-N upper triangular matrix . On output, the updated matrix .

real(kind=real64), intent(inout), dimension(:) :: u

On input, the N-element vector . On output, the rotation sines used to transform to .

private pure subroutine cholesky_rank1_downdate_cmplx(r, u)

Computes the rank 1 downdate to a Cholesky factored matrix such that . This operation only works if the new matrix is positive definite.

Arguments

Type IntentOptional Attributes Name
complex(kind=real64), intent(inout), dimension(:,:) :: r

On input, the N-by-N upper triangular matrix . On output, the updated matrix .

complex(kind=real64), intent(inout), dimension(:) :: u

On input, the N-element vector . On output, the rotation sines used to transform to .